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ABSTRACT 

We find that the orbital period (2.4 hours), eccentricity (0.09), dipole magnetic field 
strength (6.9 x 10 9 Gauss) and spin period (22 ms) of the new highly relativistic double 
neutron star system PSR J0737-3039 can all be consistently explained if this system 
originated from a close helium star plus neutron star binary (HeS-NS) in which at the 
onset of the evolution the helium star had a mass in the range 4.0 to 6.5 Mq and 
an orbital period in the range 0.1 to 0.2 days. Such systems are the post-Common- 
Envelope remnants of wide Be/X-ray binaries (orbital period ~ 100 to 1000 days) 
which consist of a normal hydrogen-rich star with a mass in the range 10 - 20 Mq 
and a neutron star. The close HeS-NS progenitor system went through a phase of mass 
transfer by Roche-lobe overflow at a high rate lasting a few times 10 4 years; assuming 
Eddington-limited disk accretion onto the neutron star this star was spun up to its 
present rapid spin rate. At the moment of the second supernova explosion the He star 
had a mass in the range 2.3 to 3.3 Mq and in order to obtain the present orbital 
parameters of PSR J0737-3039 a kick velocity in the range 70 - 230 kms" 1 must have 
been imparted to the second neutron star at its birth. 
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1 INTRODUCTION 

A new double neutron star (DNS) system, J0737-3039, was 
recently discovered (Burgay et al. 2003). This is the fifth 
galactic DNS; all share the same characteristics of having 
eccentric orbits (e > 0.08) and short pulse periods (22 - 104 
ms) and low period derivatives (~ 10 _18 ss _1 ). The two 
latter parameters imply a weak surface magnetic field (~ 
10 10 G) which indicates that the pulsar has been recycled 
(Smarr & Blandford 1976; Alpar et al. 1982; Radhakrishnan 
& Srinivasan 1982), and that the observed pulsar in each 
system is the older neutron star and has gone through a 
mass-transfer phase (Srinivasan & van den Heuvel 1982). 

Dewi & Pols (2003) found that the wide-orbit DNSs 
J15 18+4904 and J 181 1-1736 are produced from helium 
star-neutron star binaries that avoided Roche-lobe overflow 
(RLOF); while the close-orbit B1913+16 and B1534+12 can 
be formed either with or without a mass-transfer phase 
in the helium star-neutron star binary stage of evolution. 
Among the five galactic DNSs, J0737-3039 has the short- 
est spin period (22 ms), shortest orbital period (0"?102), and 
lowest eccentricity (0.088). In this paper we examine the 



email: jasinta@science.uva.nl 



possible formation history of J0737-3039. We will show that 
unlike the other DNSs, J0737-3039 can only be formed from 
helium star-neutron star binaries which have gone through 
RLOF. 



2 THE FORMATION OF J0737-3039 
2.1 Possible progenitors 

We examine the formation history of the system is terms 
of the so-called standard scenario (Bhattacharya & van den 
Heuvel 1991), in which DNSs were formed from Be/X-ray 
binaries. When the Be star leaves the main sequence and 
becomes a giant, a runaway mass transfer takes place due 
to the large mass ratio of the system. The neutron star is en- 
gulfed by the giant's envelope, and a common-envelope (CE) 
phase ensues in which the neutron star spirals down to the 
helium core of the giant, ejecting its hydrogen-rich envelope 
in this process. The outcome of this spiral-in is a very close 
binary consisting of a helium star and a neutron star (HeS- 
NS). For companions less massive than about 5 Mq, this 
mass transfer can still be stabilized such that the formation 
of a CE can be avoided (e.g. King & Ritter 1999; Podsiad- 
lowski & Rappaport 2000; Tauris, van den Heuvel, Savonije 
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Figure 1. The helium star masses and the orbital periods after 
a mass-transfer phase in HeS-NS binaries, taken from Dewi et al. 
(2002) and Dewi & Pols (2003). The open star symbols mark the 
systems which go through a CE phase due to the development 
of convective helium envelope. Thin lines connect the remnants 
of systems with the same initial helium star masses, M\. The 
shaded area marks the region where the helium stars never fill 
their Roche lobes. The maximum and minimum pre-SN orbital 
period of J0737-3039 are represented by the upper and lower thick 
lines, derived from the separation of the system at the time it was 
born, i.e. a characteristic spin-down time ago. 

2000). However, for more massive companions such as the 
Be stars considered here the formation of a CE cannot be 
avoided. 

The helium star evolves further through core-helium 
burning and further burning stages. During this evolution 
it transfers mass to the neutron star by RLOF and also 
loses mass by stellar wind. Fig. depicts the parameters of 
these HeS-NS binaries at the end of this RLOF and wind 
mass loss, resulting from detailed evolutionary calculations 
(Dewi et al. 2002; Dewi & Pols 2003). The number indi- 
cated at each curve represents the initial mass of the helium 
star at the beginning of its evolution. The stars indicate the 
outcome of the evolution of systems with a range of initial 
orbital periods (note that the lines connecting the stars are 
not evolutionary tracks). Also presented in the figure are 
the maximum and minimum possible orbital periods of the 
progenitor of J0737-3039 just prior to the supernova (SN) 
explosion. These lines are derived with the assumption that 
the pre-SN (circular) orbital separation, at, lies between the 
periastron and apastron of the post-SN separation, i.e. 

Of (1 - et) < at < at (1 + e f ) (1) 

(Flannery & van den Heuvel 1975), where at and et are the 
post-SN separation and eccentricity. Here at and et are the 
parameters at the time when it was born, which is assumed 
to be a characteristic time ago, r = P/2P. These parameters 
were derived by integrating the effects of gravitational-wave 
radiation on the observed (a, e) over r. 

Dewi & Pols (2003) found that systems represented with 
the open star symbols undergo a common-envelope (CE) 
phase due to the development of a helium convective enve- 
lope. The consequence of this is that, after the CE phase 



these systems will have a lower mass (Mn e ,t ~ 1-4 Mq) and 
much shorter orbital period (Pt ~ O^Ol) than their original 
parameters shown in Fig. Hence, only the helium stars 
represented by solid star symbols (i.e. with initial masses in 
the range of 4 - 6.5 Mq) are possible progenitors of J0737- 
3039. This range of mass requires an initial progenitor main- 
sequence star of 14 - 20 Mq. 

The region of the possible pre-SN orbital period P lim- 
ited by the two thick lines in Fig.^does not cross the shaded 
area where the helium stars never fill their Roche lobes. This 
implies that J0737-3039 can only have been formed through 
Roche-lobe mass transfer from the helium star progenitor - 
as opposed to the other four galactic DNSs which could be 
formed either via RLOF or without a mass-transfer phase 
(Dewi & Pols 2003). 

Fig. depicts an enlargement of the region of allowed 
final HeS-NS systems from Fig.^in an a t vs. MHe,t diagram 
(the stars are the same as the allowed systems in Fig.0. For 
a given combination of pre-SN helium star mass Mrie.t an d 
orbital separation at we calculated the kick velocity required 
to produce the post-SN orbital period and eccentricity. Fig.|2] 
shows the ranges of the required kick velocities V\ limited 
by the curves of V k = 70, 125, 175 and 230 kms~\ The 
region where the systems undergo a CE phase and evolve 
to extremely short orbital periods is to the left of the dash- 
dotted line in Fig. [5] Assuming that only systems which do 
not go through the CE phase are the possible progenitors 
of J0737-3039 (i.e. the systems on the right-hand side of 
the line), we find that to explain the orbital parameters of 
J0737-3039 when it was born (i.e. a — 1.688 Rq and e = 
0.138; assuming that the pulsar's age equals its characteristic 
age), a minimum kick velocity of 70 km s" 1 is required at 
the birth of the pulsar's companion. 

The minimum separation, below which mass-transfer 
phase from the helium star to a neutron star is dynamically 
unstable resulting in a merger of the two stars (Dewi et al. 
2002), is marked by the dotted line in Fig. 0] This line gives 
a constraint on the maximum kick velocity, i.e. 230 kms -1 . 



2.2 Magnetic field and spin history of the pulsar 

All five double neutron star systems (DNSs) in eccentric 
orbits in the galactic disk share the characteristic of a com- 
bination of a weak dipole magnetic field (in the range 10 9 to 
2 x 10 10 G) and rapid spin, indicating that we are dealing 
here with recycled pulsars that have been spun up during an 
accretion phase (Smarr & Blandford 1976; Alpar et al. 1982; 
Radhakrishnan & Srinivasan 1982; Srinivasan & van den 
Heuvel 1982). During such a phase with a non-degenerate 
donor companion the orbit will have been tidally circular- 
ized. Therefore the present eccentricity of these systems can 
only have been induced in the second supernova explosion 
when the companion collapsed to a neutron star (Srinivasan 
& van den Heuvel 1982) . The reduction of the magnetic field 
strength of the pulsar to a value around 10 9 to 10 10 G is 
thought to somehow be the result of the accretion, as there 
is an observed correlation between the amount of matter ac- 
creted and the field strength (Taam & van den Heuvel 1986). 
A simple relation for this case was derived semi-empirically, 
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Figure 2. The possible pre-SN separation of J0737-3039 is lim- 
ited between the two horizontal lines. To the left of the dash- 
dotted line systems undergo a CE phase and may evolve to ex- 
tremely short orbital periods. The dotted line shows the minimum 
separation allowed for the a helium star to undergo a stable mass- 
transfer phase. To the right of this line runaway mass transfer oc- 
curs, probably leading to coalescence. The solid circle marks the 
pre-SN mass and separation in case of a symmetric SN explosion 
at the birth of the pulsar's companion. The thin lines give the 
magnitude of the kick velocity (the corresponding kick velocities 
in kms -1 are presented at the bottom of the figure). 



by Shibazaki et al. (1989): 
Bo 



B(t) 



1 + AM/tub 



(2) 



where B(t) and -Bo are the present and initial magnetic field 
strengths, respectively, AM is the amount of matter ac- 
creted and me is a constant to be derived from the observa- 
tions. Francischelli et al. (2002) found m B = 1.25 x 10~ 5 Mq 
to agree well with the observed data for Low Mass X-ray Bi- 
naries. Adopting this value and a duration of the Roche-lobe 
overflow phase of the helium star of 2 x 10 4 yrs (this is about 
the average value for helium stars in the mass range 4.0 to 
6.3 M Q (Dewi et al. 2002; Dewi & Pols 2003)), with accre- 
tion at an Eddington rate for helium (3 x 10 -8 Mq yr _1 ) one 
finds that B(t) can have decreased by a factor of about 50 
during this phase alone. Furthermore, a 15 Mq Be star re- 
sides 10 7 yrs on the main-sequence. Be/X-ray binaries tend 
to be transients (Rappaport & van den Heuvel 1982; van den 
Heuvel & Rappaport 1987; Zioikowski 2002) and are "on" as 
an X-ray source only part of the time, during which they of- 
ten reach a luminosity close to the Eddington limit. Assum- 
ing them to be "on" only one per cent of the time and then to 
reach the Eddington limit, the neutron star will during the 
Be/X-ray binary phase have accreted about 1.5 x 10 -3 Mq. 
With eq. @ this is already sufficient to reduce its mag- 
netic field by a factor of about 100. So, adding the amount 
of matter accreted during the Be/X-ray binary phase and 
the RLOF of the helium star, one finds with eq. that 
these phases together provide enough accretion to reduce 
the magnetic field by a factor of about 200. 



Magnetic fields of young pulsars on average typically 
have a dipole strength a few times 10 12 G, and after re- 
duction by this factor will typically be around 10 10 G, as 
observed in the five DNSs. As to the spin-up of the "old" 
neutron star: during the Roche-lobe overflow from the he- 
lium star accretion to the neutron star will have proceeded 
through a disk. The amount of angular momentum accreted 
during this phase can be calculated as follows. The mass 
transfer rate during Roche-lobe overflow is highly super- 
Eddington, in the range 1 — 2 x 10~ 4 Mq yr _1 for helium 
stars in the mass range 4.0 - 6.3 Mq (Dewi et al. 2002; 
Dewi & Pols 2003). Dewi et al. (2002) and Dewi & Pols 
(2003) in their calculations of the evolution of the systems 
during this phase assumed the accretion rate onto the neu- 
tron star to be the Eddington rate (3 x 10~ 8 Mq yr _1 ) and 
the super-Eddington fraction of the transferred matter to 
be blown out of the system by the radiation pressure, car- 
rying with it the specific angular momentum equal to the 
specific orbital angular momentum of the neutron star. Fol- 
lowing this assumption of accretion at the Eddington rate 
- as radiation pressure will prevent accretion at a higher 
rate - and using the dipole strength of the magnetic field of 
J0737-3039 of 6.9 x 10 9 G, one can calulate the value of the 
Alfven-radius R R , where the matter becomes coupled to the 
magnetic field and starts to co-rotate with the neutron star 
magnetosphere. 7? a is given by (e.g. see Bhattacharya & van 
den Heuvel 1991)): 



{BsRlf 7 



(M/Me,h) 2 / 7 (M/M Q )!/7 



.32.9 km 



(3) 



where Bg is the dipole field strength in units of 10 9 G, Re 
is the neutron star radius in units of 10 6 cm, M and Me,h 
are the accretion rate and the Eddington accretion rate for 
hydrogen, respectively and M is the mass of the neutron 
star. With the dipole field strength of 6.9 x 10 9 G of PSR 
J0737-3039, the Eddington rate for helium accretion M = 
M E ,He = 3 x 10~ 8 Mq yr" 1 and a M = 1.3 Mq one obtains 
i? a = 78.4 km. 

The rate at which angular momentum is accreted by 
the neutron star is then: 



J = MQ,*(r = R 3 )Kl 



(4) 



where fik(>") is the keplerian angular velocity of matter 
in the disk at distance r from the center of the neutron 
star. With M = M E ,Ho, M = 1.3 Mq one finds j 
1 , 1 x 10 _1 Mq km 2 s _1 yr _1 . The present spin angular mo- 
mentum of the neutron star is 

J = M(2ty/P) R 2 r'l (5) 

where R and P are the radius and spin period of the neutron 
star and r g is its radius of gyration, respectively. The value of 
rt is typically about 0.075 for a polytropic index n — 3 (e.g. 
Lattimer & Prakash 2001). With this value and R = 10 km 
and P = 22 ms one obtains J = 2.8 x 10 3 M Q km 2 s~ 1 , 
which together with the above mentioned rate of accretion 
of angular momentum results in a required timescale for the 
spin-up of about 2.5 x 10 4 years. 

The duration of the phase of Roche-lobe overflow in 
the helium star plus neutron star binaries with helium star 
masses between 4.0 and 6.3 Mq ranges from 2.4 to 0.7 x 
10 4 years, which is of the same order as the value of 2.5 x 
10 4 years required for the spin-up. It thus appears that the 
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magnetic field strength and spin rate of PSR J0737-3039 fit 
quite well with evolutionary model presented here in which 
the immediate progenitor system was a very close helium 
star plus neutron star binary which evolved through a phase 
of Roche lobe overflow with a duration of order 10 4 years. 
As mentioned above, the progenitor of this system was so 
close that it cannot have avoided the phase of mass transfer 
by Roche-lobe overflow from the helium star (Dewi & Pols 
2003), so here this phase must absolutely have taken place in 
the history of PSR J0737-3039. This may be the reason why 
the pulsar has the shortest spin period of all. The progenitors 
of the other systems may, in principle, have avoided RLOF, 
although it can certainly not be excluded that also in (some 
of) them Roche-lobe overflow from the helium star may have 
taken place. 



3 DISCUSSIONS AND CONCLUSIONS 

We have examined the formation of J0737-3039. We find 
that this system must have gone through a mass-transfer 
phase from a helium star with the initial mass in the range 
of 4 - 6.3 Mq . We also find that an asymmetric SN explosion 
with a kick velocity of 70 - 230 km s _1 is necessary to explain 
the observed orbital parameters. However, these constraints 
may still be subject to some uncertainties. Whether the neu- 
tron stars in the systems that go through a CE phase (the 
open stars in Fig. complete the spiral-in in the envelope 
of the helium star depends on the competition between the 
spiral-in timescale and the time left before the explosion of 
the helium star (Dewi & Pols 2003). If the SN explosion in 
systems with open star symbols in Fig. Q occurs before the 
neutron star completes the spiral-in process, a lower kick ve- 
locity or even a symmetric SN explosion (shown as the solid 
circle) would then be possible. Furthermore, the constraints 
on the maximum pre-SN mass and the maximum kick veloc- 
ity depend somewhat on the assumed wind mass-loss rate 
and the loss of angular momentum during the RLOF. As the 
stability of a mass-transfer phase depends on the mass ratio 
of the system and the period at the onset of RLOF, applying 
a higher wind mass-loss rate results in a lower mass helium 
star in wider orbit, and hence, more stable mass transfer 
than adopting wind mass loss with a lower rate. 
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